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Abstract 



In quantum field theory, heavy particles with mass M much greater than the tem- 
perature T give not only effects suppressed by the Boltzmann factor e~ M l T , but 
also effects suppressed by powers of T/M. We show that power-suppressed terms in 
equilibrium observables arise from effective interactions among light particles due to 
virtual heavy particles. We study a model introduced by Matsumoto and Yoshimura 
in which heavy bosons interact only through a term that allows pair annihilation into 
light particles. We construct an effective Lagrangian for the light field by integrating 
out the heavy field, and use it to calculate the leading power-suppressed terms in the 
energy density. The thermal average of the Hamiltonian density for the heavy field 
includes a term proportional to T 6 /M 2 , but we show that this term can be eliminated 
by a field redefinition and therefore cannot have any physical significance. 
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I. INTRODUCTION 



One might naively expect the effects of heavy particles with mass M much greater than 
the temperature T to be suppressed by the Boltzmann factor e~ M l T . However, in a quantum 
field theory, there are additional effects that are suppressed only by powers of T/M. In a 
recent series of papers |T|-^|, Matsumoto and Yoshimura have argued that there are power- 
suppressed terms in the number density of heavy particles. If there were such terms, they 
could greatly exceed the contribution from the conventional Boltzmann-suppressed terms 
when T <C M. This would have important implications for cosmology, because it would 
imply that the relic abundance of weakly-interacting massive particles is much larger than 
the conventional predictions based on the Boltzmann equation. Present bounds on the 
energy density of the universe would then imply significantly tighter constraints on the 
properties of the heavy particles that may constitute the cold dark matter. 

Matsumoto and Yoshimura have studied the power-suppressed effects in a simple model 
with two species of scalar particles, one heavy and one light. The only interaction of the 
heavy particle is one that allows pair-annihilation into light particles. In the first two papers 
in the series Matsumoto and Yoshimura used the influence functional method and a 
Hartree approximation to derive a quantum kinetic equation for the momentum distribution 
of heavy particles. Their equation includes off-shell effects associated with the thermal width 
of the heavy particles. They found that the leading power-suppressed contribution to the 
heavy-particle number density was proportional to T 7 / 2 / 'M 1 / 2 . In their third paper M, 
Matsumoto and Yoshimura studied the equilibrium number density of heavy particles and 
found that the leading power-suppressed term was actually proportional to T 6 /M 3 . In 
their fourth paper [Q], they derived a new quantum kinetic equation that reproduces the 
equilibrium result of Ref . Q . 

Singh and Srednicki M have criticized Matsumoto and Yoshimura's conclusion that there 
are power-suppressed contributions to the number density of heavy particles. They argued 
that the quantum kinetic equation of Matsumoto and Yoshimura does not properly account 
for the interaction energy between the heavy particles and the thermal bath of light particles. 
They also noted that the power-suppressed terms found by Matsumoto and Yoshimura are 
actually contributions to the number density of virtual heavy particles. They argued that 
it is the number density of on-shell heavy particles that is relevant to the relic abundance, 
and this will have the usual Boltzmann suppression. 

Srednicki recently argued that Matsumoto and Yoshimura's conclusion is the result 
of an inappropriate definition of the number density ||. Srednicki considered a similar 
model in which the real-valued heavy field is replaced by a complex-valued field, so that 
the heavy bosons have a conserved charge. He showed that there was a definition of the 
heavy-particle number density that had Boltzmann suppression to all orders in perturbation 
theory. Srednicki also suggested that the power-suppressed contributions to the heavy- 
particle energy density should have a simple interpretation in the effective field theory for 
the light particles obtained by integrating out the heavy particles. They are contributions 
to the energy density of light particles coming from nonrenormalizable effective interactions 
between the light particles. 

The approach followed by Matsumoto and Yoshimura has been to integrate out the light 
field to get an quantum kinetic equation for the heavy particles. The power-suppressed terms 
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are then understood as arising from the thermal width acquired by the heavy particle when 
the light field is integrated out. The heavy particle no longer has a sharp energy-momentum 
relation, but is a resonance. The power-suppressed terms come from the tail of the spectral 
function of the resonance, where the energy and momentum of the heavy particle are both 
small compared to M. 

The philosophy of effective field theories suggests using the diametrically opposite 
strategy. Physics involving energies and momenta small compared to M can be understood 
most simply by integrating out the heavy field. A heavy particle whose energy and mo- 
mentum is small compared to M is off its mass-shell by an amount of order M. By the 
uncertainty principle, it can remain in this highly virtual state only for a time of order 1/M. 
Light fields can propagate only over distances of order 1/M in this short time. Thus the 
effects of the highly virtual heavy particle on light fields with momenta much smaller than 
M can be taken into account through local interactions among the light fields. In other 
words, the light fields can be described by a local effective field theory. 

The effective field theory approach was used by Kong and Ravndal |7j] to compute the 
leading power-suppressed terms in the energy density for QED at temperature T <C m e , 
where m e is the electron mass. They first integrated out the electron field to get a low- 
energy effective Lagrangian for photons that includes the Euler-Heisenberg term. They then 
computed the energy density for this effective field theory at temperature T and found that 
the leading term is proportional to a 2 T s /m^. This term can be identified as a contribu- 
tion to the photon energy density coming from the Euler-Heisenberg term in the effective 
Hamiltonian for photons. 

In this paper, we use similar effective field theory methods to study the power-suppressed 
thermal effects in the model considered by Matsumoto and Yoshimura. We show that the 
power-suppressed terms in the energy density can indeed be interpreted as contributions to 
the energy density of light particles from nonrenormalizable effective interactions. We also 
show that the term in the energy density from which Matsumoto and Yoshimura extracted 
the heavy-particle number density can be eliminated by a field redefinition and therefore 
can not have any physical significance. 

This paper is organized as follows. We introduce the pair-annihilation model of Mat- 
sumoto and Yoshimura in section II and summarize their results on contributions to the 
energy density that are suppressed by powers of T/M. In section III, we construct an effec- 
tive Lagrangian for the light field by integrating out the heavy field. We compute the leading 
power-suppressed terms in the energy density by differentiating the free energy density for 
the effective theory in equilibrium at temperature T. In section IV, we use a field redefinition 
to construct an effective Hamiltonian for the light field. We show that the thermal average 
of the effective Hamiltonian density reproduces the leading power-suppressed terms in the 
energy density. We summarize our results in section V. 

II. BOSONIC MODEL WITH PAIR ANNIHILATION 

The model studied by Matsumoto and Yoshimura contains two species of spin-zero 
particles: a heavy particle of mass M described by the field (p and a massless particle 
described by the field \. The Lagrangian density is 
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c = c x + \d^d^ - ImV - i W, (i) 

where C x is the Lagrangian for the light field: 

£x = \d&Px - ^A xX 4 . (2) 

We have suppressed the counterterms needed to remove ultraviolet divergences. The sym- 
metry (f — > — y> guarantees that the heavy particle is stable at zero temperature. The heavy 
particles can be created or annihilated in pairs via the <p 2 x 2 interaction. Renormalizability 
also requires a ip A self-interaction, but we assume that its coefficient is much smaller than 
A, so it can be neglected. This is a great simplification, because the Lagrangian is then 
quadratic in the heavy field (p. The x 4 term in (P is necessary for thermalization of the 
light field, but since we are primarily interested in the effects of the heavy field, we will carry 
out explicit calculations only to zeroth order in A x . 

The energy density is the ensemble average of the Hamiltonion density: p = (Ti). Mat- 
sumoto and Yoshimura divide the Hamiltonian density into three terms: Tt = T~t x +TCtp+TCi n t, 
where 

^x = ^ 2 + ^(Vx) 2 + ^A xX 4 , (3) 

^ = \^ + \{V^ + \m^\ (4) 

H mt = \\^\ 2 . (5) 

We have suppressed the counterterms required to renormalize the composite operators so 
that their expectation values vanish at zero temperature. Matsumoto and Yoshimura in- 
terpreted the corresponding three terms in p = p x + p v + p int as the energy density of the 
"thermal environment", the energy density of the "system" consisting of heavy particles, 
and the interaction energy density, respectively. 

At zeroth order in A, the energy density of the heavy particles is that of an ideal non- 
relativistic gas, 

p v = M{MT/2itf' 2 e- M / T , (6) 

which exhibits the usual Boltzman suppression. At second order in A, there are terms that 
are suppressed only by powers of T/M. Matsumoto and Yoshimura calculated the leading 
power-suppressed contributions for each of the 3 terms in the energy density @: 

1 T 6 

Sf>y = A 2 -^> (7) 

Hx 69120 M 2 v ; 

1 T 6 

SPv = 69120 A V' (8) 

5p ' mt = ~64^00 A 2 W (9) 

The terms proportional to A 2 T 6 terms cancel between 5p x and 5p v and, so the leading 
power-suppressed term of order A 2 in the total energy density is proportional to A 2 T 8 . 
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Matsumoto and Yoshimura noted this cancellation, but nevertheless interpreted 5p v in (§) 
as a contribution to the energy density of heavy particles. Taking the heavy particles to 
be nonrelativistic with energy equal to M, they identified 5p v /M as a contribution to the 
heavy-particle number density proportional to T 6 /M 3 . 

Singh and Srednicki Q have argued that the separation of the energy density into three 
terms corresponding to the system, the environment, and interactions is reasonable if and 
only if | p int | <C p^. If this condition is not satisfied, the coupling between the system and 
the environment is effectively strong and they cannot be clearly separated. Note that this 
condition is satisfied by the power-suppressed terms (§) and (|S|) if T <C M. 



III. EFFECTIVE LAGRANGIAN 



In this section, we construct a low-energy effective Lagrangian for the light field \. This 
effective Lagrangian reproduces the zero-temperature Green functions at momentum scales 
much less than M. The effects of the heavy fields are reproduced by nonrenormalizable 
interactions with coefficients that are suppressed by powers of 1/M 2 . 

An effective Lagrangian for the light particles can be constructed by using functional 
methods to integrate out the heavy field. This method is particularly convenient for the 
model of Matsumoto and Yoshimura, because the Lagrangian is quadratic in the field <p. 
The effective action for the light field can be defined by a functional integral over the heavy 
field: 

exp(iS c g[x]) = J Vipexp(i J d 4 x C). (10) 

This effective action shouldn't be confused with the 1PI effective action that generates one- 
particle-irreducible Green functions. Since £ is quadratic in <p, the functional integral can 
be evaluated explicitly: 

S cS [ X ]= Jd A xC x + l -lndet^-d 2 -M 2 -^x 2 + ^y (11) 

The effective action can be expanded in powers of the coupling constant A: 

S cS [ X ]= d A x£ x + -\ndet(-d 2 -M 2 + l e) + J2S$ix}- (12) 

n=l 

The term of n'th order in A is 



3^ \x\ = "^tr [(-d 2 — M 2 + *e)-Vj " . (13) 



The lndet(-<9 2 - M 2 ) term in flU) 

can be discarded, because it is just a ^-independent 
constant. The only contribution from S^g is the local functional / d 4 xx 2 with a divergent 
coefficient. It can be absorbed into the mass counterterm for the light field. The terms S^g 
for n > 2 are nonlocal functionals of the x field. 
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Since we are interested in light fields with characteristic momenta much smaller than M, 
we can use the derivative expansion to express S^g as an infinite series of local functionals. 
The derivative expansion is illustrated in Appendix A by expanding S^ s to all orders. The 
lowest derivative term from is a x 4 term. It has a divergent coefficient and can be 
absorbed into the counterterm for the x 4 term in the Lagrangian. The remaining terms 
have finite coefficients suppressed by powers of 1/M 2 . They represent nonrenormalizable 
interactions among the light fields induced by virtual heavy particles. The effective action 
can now be expressed as the integral of an effective Lagrangian: S^rtx] = / d A xC e s, where 

£„ = c — y 2 <9V + — y 2 (d 2 ) 2 y 2 — y 6 + (U) 

cS x 96(47r) 2 M 2 * X + 960(4tt) 2 M 4X 1 ' X 96(4tt) 2 M 2X + "" { } 

We have suppressed the counterterms required to remove ultraviolet divergences, and we 
have kept all terms proportional to A m (l/M 2 ) n with m + n < 4. If we consider an observable 
involving a single momentum scale p <C M, terms in £ c g proportional to A m (l/M 2 ) n will 
give effects suppressed by X m (p/M) 2n . The x 6 term in (|TJ]) will therefore be comparable in 
importance to the x 2 {d 2 ) 2 X 2 term if A ~ (p/M) 2 . 

A similar strategy can be used at nonzero temperature T to compute power-suppressed 
contributions to equilibrium observables if T -C M. Such observables can be expressed as 
Euclidean functional integrals over the fields ip and x with periodic boundary conditions 
in the Euclidean time direction. To calculate the power suppressed terms, we would first 
integrate over ip, then expand in powers of A, then carry out the derivative expansion in 
powers of 1/M 2 , and finally integrate over x- The first three steps reduce the problem to a 
calculation in an effective theory for x, and the final step of integrating over x corresponds 
to computing the thermal average in that effective theory at temperature T. This strategy 
will reproduce all the terms that are suppressed by powers of T/M, but it will not give any 
Boltzmann-suppressed terms, because the expansion in powers of 1/M 2 eliminates terms 
with an essential singularity at 1/M = 0. 

The simplest way to compute the power-suppressed terms in the energy density p is to 
calculate the corresponding terms in the free energy density T and then differentiate. The 
power-suppressed terms in T can be obtained simply by computing the free energy density 
at temperature T for the theory defined by the effective Lagrangian (|l~4"D. The leading terms 
are given by vacuum diagrams whose only vertex is one of the power-suppressed interactions 
in ([14]), and they can be written 
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bT 96(4tt) 2 M 2 ^ 2<92 * 2 > frcc 960(47r) 2 M 4 

+ 96(4^) 2 M 2 (xVe + -- - ' (15) 

The angular brackets (...)f ree denote the thermal average in the free field theory. These 
thermal averages are expressed as Matsubara sum-integrals in Appendix B. The first term 
on the right hand side of (JI5D is zero. The remaining two terms give 



1 / 16 T 4 25 T 2 \ 
^ = 1024 225 A V ~ 48^ A V + ' ' "J ^ (16) 
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where J-f ree = — (7r 2 /90)T 4 is the free energy of a gas of free massless bosons. By dimensional 
analysis, the x 2 9 2 X 2 term in fll5D would have given a term proportional to A 2 T 6 /M 2 . The 
absence of such a term is related to the cancellation of the A 2 T 6 /M 2 terms in the energy 
density noted by Matsumoto and Yoshimura. 

Once the free energy density is known, we can derive the energy density by differentia- 
tion: p = -T 2 ^(^/T). The leading power-suppressed terms in the energy density are 

1 /112 2 T 4 _ 125 3 T 2 
P ~ 1024 V675 M 4 ~ 144vr 4 M 2 

where pf ree = (7r 2 /30)T 4 is the energy density of a free gas of massless bosons. Srednicki 
calculated the A 2 (T/M) 4 term for a similar model with a complex- valued heavy field by 
direct calculation in the full theory ||. It is worth noting that the A 3 (T/M) 2 term is equally 
important if A ~ (T/M) 2 . At very low temperature, the power-suppressed terms in fll7j) 
dominate over the leading Boltzman-suppressed term (|g) in the energy density of the heavy 
particles, but they represent small corrections to the energy density of the light particles. 




IV. EFFECTIVE HAMILTON! AN 

Srednicki [|| argued that the power suppressed terms in the energy density should be 
interpreted as contributions to the energy density of the light field from nonrenormalizable 
effective interactions. In order to verify this explicitly, we construct a low-energy effective 
Hamiltonian density for the light field and compute its thermal average in the effective 
theory at temperature T. 

If the Lagrangian density depends only on the field x an d its first derivatives, the 
standard Noether prescription for constructing the Hamiltonian density is TC = x{dC/dx) — 
C. This prescription can not be applied to the effective Lagrangian fll4]) because, even 
after using integration by parts to reduce the number of derivatives acting on any single 
field, it still depends on the second derivatives of x- While the Noether prescription for 
the Hamiltonian density can be generalized to higher derivative Lagrangians, it is rather 
cumbersome. 

A simpler approach is to first construct a different effective Lagrangian C' eS that depends 
only on x and its first derivatives, and then apply the Noether prescription to it. The effec- 
tive Lagrangian £ eff in (O) is the unique effective Lagrangian that reproduces the off-shell 
Green's functions of the full theory at low momenta, but there are infinitely many effective 
Lagrangians that reproduce all the physical observables at low momenta. They include all 
effective Lagrangians that can be obtained from (|14]) by a field redefinition. In quantum field 
theory, we always have the freedom to redefine the field, because physical quantities, such 
as S-matrix elements, are invariant under field redefinitions. In renormalizable field theo- 
ries, nontrivial field redefinitions are usually not considered, because they make the theory 
superficially nonrenormalizable. However, effective theories are already nonrenormalizable, 
so nontrivial field redefinitions don't introduce any additional complications. In fact they 
can be used to simplify the effective Lagrangian by removing terms that don't contribute to 
physical quantities. For example, by introducing the field redefinition x - * X + G(x) hito 
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the kinetic term d^x^Xi we generate additional terms that can be used to cancel any terms 
of the form G(x)d 2 x- 

We use the following field redefinition to simplify the effective Lagrangian (|14"D: 

v > v y 3 + dV + ... . (18) 

A A 72(4vr) 2 M 2A 720(4tt) 2 M 4 A 1 J 

Expanding out the derivatives and rearranging them by using integration by parts, our 
effective Lagrangian reduces to 

A 2 A 3 
C ' cS = Cx + 2A^fM^ x9 " x)2 ~ 96(4tt) 2 M 2x6 + " ' ' (19) 

Again we have kept only those terms with a total of up to 4 powers of A and 1/M 2 . This new 
effective Lagrangian will not reproduce the low-momentum Green functions of the original 
theory, but it will reproduce all physical observables involving low momenta. 

An alternative way to derive the effective Lagrangian (|19|) starting from the original 
Lagrangian ([!]) is by matching physical quantities computed in both theories. We would 
begin by writing down the most general effective Lagrangian consistent with the symmetry 

X -> -x- 

C' cS = C X + A x 2 d^X + B X 6 + C^xd^X? + D x 2 d 2 X d 2 X + E x d 2 X d^X + ■■■ ■ (20) 
We would then write down the most general field redefinition consistent with the symmetry: 
X -> X + ad 2 X + b X 3 + c X d ilX d ll X + d X 2 d 2 X + ■■■ ■ (21) 



Inserting this field redefinition into (20) and expanding it out, we would find that the 
coeffcients b, c, and d could be used to set A = D = E = 0. To determine the remaining 
coefficients, such as B and C, we would exploit the fact that physical quantities are invariant 
under field redefinitions. We would compute T-matrix elements involving light particles with 
momenta p -C M in the full theory using the original Lagrangian (|1|) and in the effective 
theory using the effective Lagrangian (20). By matching these T-matrix elements, we would 
deduce that B and C have the values given in (]TTj|) . 

Having constructed the new effective Lagrangian C' eS in (|T9|) that depends only upon x 
and its first derivatives, we can use the standard Noether prescription to deduce the effective 
Hamiltonian. The effective Hamiltonian density reads TC e s = 7~i x + TC pQ w, where H x is given 
in @ and TC pow includes all the higher dimension operators: 

«~ = mSfM' (a ^ ( 3 * 2 + (v * )2 ) + n^w^ + - • (22) 

We can now calculate the power-suppressed terms in the energy density by taking the 
thermal average (H e s) at temperature T for the effective theory defined by dl9|). The energy 
density can be written as p = (Ti x ) + (Hpow)- 111 (T^pow), the leading power-suppressed terms 
Sppo W are simply the thermal averages in a free field theory of the operators in (|22]). In (7i x ), 
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the leading power-suppressed terms Sp x come from treating the interactions in fll9|) as first- 
order perturbations. The calculations are described in more detail in Appendix B, and the 
results are 



a 1 fi6 x2 r 4 

Px 1024 I 135 M 4 



25 T 2 



5p 



pow 



1024 I 675 M 4 



24vr 4 
25 



M2 J Pfree, 

A 3 * 



144tt 4 " M 2 "i Pfree ' 



(23) 
(24) 



where pf ree = (7r 2 /30)T 4 . The sum of ( p3[) and (^) reproduces our previous result (|lj 

differs from the power-suppressed terms in (H x ) in the full theory, 



Note that Sp x in ([ 

which include the term (0) suppressed by A 2 T 2 /M 2 . 

If we had omitted the A 2 x 3 term in the field redefinition (|l8l), the term proportional to 
X 2 d 2 x 2 m the effective Lagrangian (|T4| ) would not have been eliminated. There would then 
have been an additional term in the effective Hamiltonian proportional to x 2 {x 2 + i^x) 2 )- 
Its thermal average reproduces the term 5p v in (|^) calculated by Matsumoto and Yoshimura. 
However the x 2 9 2 X 2 interaction term in the effective Lagrangian gives an additional term in 
(H x ) that reproduces Sp x in (^). Since the cancelling contributions (0) and (|8|) to the energy 
density can be eliminated by a field redefinition, neither can have any physical significance. 
In particular, 5p x /M cannot be interpreted as a contribution to the number density of heavy 
particles. 

A field redefinition was also used by Kong and Ravandal in their calculation of the 
energy density for QED at T m e . The effective Lagrangian obtained by integrating out 
the electron field is 



eff 



a 



-F, u d 2 F^ + 
607rm^ 



a 



90m 4 



{F, v F^f + 7 -(F, l/ Pn 2 



+ ... 



(25) 



The two power-suppressed terms are called the Uehling term and the Euler-Heisenberg term, 
respectively. The Uehling term can be eliminated by a field redefinition: 



-4, 



a 

M+ 3CW 



d l A^ + ... 



(26) 



It therefore cannot contribute to physical quantities. The leading power-suppressed term in 
the energy density comes from the Euler-Heisenberg interactions. 



V. CONCLUSION 



The effective- field-theory approach provides a simple way of understanding the contri- 
butions to equilibrium observables that are suppressed by powers of T/M. They arise from 
effective interactions among light particles that are induced by integrating out virtual heavy 
particles. The most economical way to compute the power-suppressed terms is to first con- 
struct a low-energy effective Lagrangian that describes the light particles at T = and then 
consider this effective theory in equilibrium at temperature T. 



9 



For the pair annihilation model introduced in Refs. [|l]-[J , we demonstrated explicitly that 
the power-suppressed terms in the energy density can be interpreted as contributions from 
the light particles. We used the field redefinition fllB|) to construct an effective Hamiltonian 
density 7i C R for the light field x> an d then verified that its thermal average reproduces the 
power-suppressed terms. The field redefinition eliminated terms suppressed by A 2 (T/M) 2 
from individual terms in (7i e s), which otherwise would have canceled only after all such 
terms had been added together. The fact that all the A 2 (T/M) 2 terms can be eliminated 
by a field redefinition indicates that individual terms of this form cannot have any physical 
significance. 

The incorrect conclusions concerning the heavy-particle number density in Ref. 0] stem 
from the authors having interpreted Ti^ in @) literally as an operator that creates only heavy 
particles and whose thermal average therefore probes the number density of those particles. 
However, the operator Ti,^ also creates light particles through loop diagrams that involve 
virtual heavy particles. Provided the momenta of the light particles are small compared to 
M, the loop diagram can be expressed as the product of a short-distance coefficient that 
depends on M and a local effective operator that creates light particles. In the definition 
of the composite operator (|]), the terms with effective operators x 2 ; X 2 ; (Vx) 2 and x 4 
are implicitly subtracted, so Ti^ creates light particles through higher dimension effective 
operators. Thus (7~C<p) receives contributions not only from heavy particles, but also from 
light particles created by these effective operators. It is these latter contributions that are 
responsible for the power-suppressed terms in the energy density. Those terms cannot be 
related to the number density of heavy particles that can participate in kinetic processes, 
since they involve only virtual heavy particles with lifetimes of order 1/M. 

The strategy of effective field theory is to integrate out heavy fields to get an effective 
theory for light fields. The construction of a quantum kinetic equation for heavy particles 
requires exactly the opposite strategy. Light fields must be integrated out to create an 
effective description of the heavy particles. Effective field theory demonstrates convincingly 
that the quantum kinetic equations derived in Ref. [3] do not describe correctly the evolution 
of the number density of heavy particles. Perhaps the insights from effective field theory 
can be used as guidance for deriving the correct quantum kinetic equations. 
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APPENDIX A: THE DERIVATIVE EXPANSION 

The derivative expansion can be used to express each of the term S^g in the effective 
action (|T2|) as an infinite series of local functionals. In this appendix, we illustrate the 
derivative expansion by applying it to the term Sjg [x]- The operator (— d 2 — M 2 + ie)~ l in 
the definition ( |13"D of corresponds to the free spin-zero propagator: 
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G(x,y) 



d 4 q 



-iq-(x-y)_ 



M 2 



IE 



(A.l) 



The definition for can be written 



S%\x] = --^J d ' x J d 4 yG(x,y)x 2 (y)G(y,x) X 2 (x). 
Inserting the integral expression ( |A.1|) for the propagators, this becomes 



(A.2) 



d 4 p 

W) 4 



ip-(x-y) 



Hp 2 



(A.3) 



where the function I(p 2 ) is given by 

d 4 q 



HP 2 



(2vr) 4 q 2 -M 2 + ie (q + p) 2 — M 2 + ie 



(A.4) 



This integral has a logarithmic ultraviolet divergence that can be isolated by adding and 
subtracting 1(0). The difference between the two integrals is convergent and can be evaluated 
using the Feynman parameter method: 



I(p 2 ) = 1(0) 



(4tt) 2 Jo 



dx In 



1 - x(l - x)p 2 /M' 



(A.5) 



Assuming the integral in (|A.3 ) is dominated by p 2 -C M 2 , we can expand the logarithm into 
a power series in p 2 and then evaluate the Feynman parameter integral to get 



Hp 2 ) = m 



E 



n\(n — 1)! / p 2 



(4tt) 2 ^ (2n + l)! V^" 2 , 



(A.6) 



The function I(p 2 ) inside the integral over p in ( |A.3|) can be replaced by I(—d 2 ) outside the 
integral. The integral over p then reduces to S A (x — y), which collapses the expression to an 
integral over a single coordinate x. Our final result for the derivative expansion is 



•\2 r \2 oo 



-l) n n\(n- 1)\ 



16(47r) 2 ^ (2n + l)!M 2r ' 



d*x X 2 (d 2 ) n x 2 



(A.7) 



The x 4 term has a divergent coefficient, but it can be absorbed into one of the counterterms 
in the original Lagrangian. All the higher derivative terms have finite coefficients. 



APPENDIX B: THERMAL SUM-INTEGRALS 

The calculations of thermodynamic quantities in this paper can be reduced to computing 
thermal averages in a free field theory. In the imaginary-time formalism, these thermal 
averages are expressed as sums over Euclidean energies and integrals over spacial momentum. 
We use the following notation for these sum-integrals: 
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The Euclidean 4-momentum is P = (p,p4 = 2irnT), where n is any integer. We also use the 
notation P 2 = p 2 + p\. 

Many of the sum-integrals required in this paper are thermal averages of local operators 
in a free field theory: 

(xW) fr ee = 4^ Q i (B.2) 



1 



3 



(xVc = 15(^ p -^J , (B.4) 
m,\0"\)' 2 h™ = t t P2Q2+ J^' Q) \ (B.5) 



P 2 Q 

((^x)(Vx) 2 ) frcc = -f p f Q PW+ p 2 P Q2 QP ' q , (B.6) 

We also need several sum-integrals that come from computing the thermal average of the 
free-field Hamiltonian density to first order in the effective interactions: 

(H X i J (fx ( W V. = -2g f Q '^^^^r 2 - 2 " 2 ' . (B.8, 

(H x i J <?* X %„ = -45 ij) 2 £ (B.9) 

(«x * / d'x X % x a" x ) bc , = ^^ (p2 + g^ 2 " 2p2) . (B.10) 



We can simplify the sum integrals by averaging over angles using (p l ) = and (p l p>) = 
p 2 S lj /3. This reduces the double sum-integrals to products of single sum-integrals. The 
single sum-integrals that are needed are 



p 



P 2 = 0, (B.ll) 

1=0, (B.12) 

p 

, i = — T 2 , (B.13) 

pP 2 12 ' v ; 



2 2 

" T 4 , (B.14) 



P 7T ^ 4 



pP 2 30' 



P 2 1, 



p (P 2 ) 2 8 ' 



T , (B.15) 



P_) " 77 T 4 . (B.16) 



p (P 2 ) 2 12 
12 



We have given only the temperature-dependent terms in the sum-integrals. The 
temperature-independent terms are ultraviolet divergent and depend on the choice of ul- 
traviolet cutoff. The most convenient cutoff is dimensional regularization of the integrals 
over the spatial momenta. With this cutoff, the temperature-independent terms vanish. 
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